On the observation of the spin resonance in superconducting CeCoIn$ 
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Recent observation of a resonance spin excitation at (1/2, 1/2, 1/2) in the superconducting state 
of CeCoIns [C. Stock et al., Phys. Rev. Lett. 100 087001 (2008)] was interpreted as an evidence for 
d x 2_ y 2 gap symmetry, by analogy with the cuprates. This is true if the resonance is a spin exciton. 
We argue that such description is undermined by the three-dimensionality of CeCoIn^. We show 
that in 3D systems the excitonic resonance only emerges at strong coupling, and is weak. We argue 
in favor of the alternative, magnon scenario, which does not require a d x 2_ y 2 gap. 



PACS numbers: 74.20.-z, 74.25. Gz, 74.72.-h 

Whether Heavy Fermion materials (HF) and high 
temperature superconductors (HTSC) are governed by 
the same physics has been repeatedly discussed in the 
past. Speculations about the similarities between the 
two classes of materials started immediately after Bed- 
norz and Muller's discovery of the high-temperature su- 
perconductivity in the cuprates [lj. To a large extent, the 
discussion at that time centered around an observation 
that in both systems spin degrees of freedom are among 
low-energy excitations, and superconductivity should be 
unconventional if spins are involved in the pairing. It 
was later realized that, although electron-electron inter- 
actions and magnetism do play indispensable roles in 
both systems, the HTSC are different from HF in several 
important aspects. Among them: (i) the metallic state 
in the cuprates emerges as a result of a chemical dop- 
ing into the antiferromagnetic (AFM) Mott insulator; (ii) 
the transport in the cuprates is very anisotropic, and the 
consensus is that the new physics is two-dimensional(2D) 
and uniquely related to CuOi planes; (iii) ARPES exper- 
iments in the optimally doped and underdoped cuprates 
do not see well-defined single particle excitations above 
T c , except in the so-called nodal regions [2]. On the 
contrary, in HF, like UPts or CeAl^, the de Haas-van 
Alphen (dHvA)-type experiments clearly indicate that 
the physics is three-dimensional (3D), and that single- 
particle excitations remain metallic and well-defined, 
although electronic characteristics dramatically change 
upon lowering the temperature below an effective " Kondo 
temperature" , Tr , at which the local f-electron states 
"mix up" with the s-p-d bands. 

The still unresolved issue is whether the difference in 
the effective dimensionality and the nature of normal 
state excitations gives rise to a different symmetry of the 
superconducting order parameter in the HF compared to 
HTSC. There is a consensus among researchers that the 
SC order parameter in HF is indeed non-s— wave, but 
there is no agreement about its symmetry (see Ref|3j for 
the summary of early results on HF ). The interest to 
this issue resurfaced recently after the discovery of the 
new class of HF with the chemical formula CeMIn^ (la- 



beled as CellS), where M = Co,Rh,Ir (Ref|3). These 
materials are much closer to the cuprates than other HF 
- they crystallize into the layered tetragonal structure 
and possess a considerable anisotropy in their electronic 
properties. The dHvA experiments have found that at 
least one FS (electron 15-band ) has apronounced quasi 
two-dimensional (Q2D) character [H, [if (see also Ref0). 
The transport properties of Cell5s also display surpris- 
ing similarities with Q2D cuprates [111 ]. 

Another remarkable feature of the Cell5 materials, 
also reminiscent of the cuprates, is the proximity between 
SC and AFM. The close interrelation between supercon- 
ductivity and antiferromagnetism has been demonstrated 
by the experiments on different alloys of Cell5s (Ref|8|). 
The transition from SC to AFM can be tuned by applied 
pressure or the magnetic field (see Refi and references 
therein), or by substituting Cd for In [101 ] . 

Such closeness between HTSC and Cell5s stimulated 
anew speculations that the superconductivity in both 
families of materials may have the same (magnetic) ori- 
gin, and the gap symmetry is the same However, un- 
like HTSC, where the d x 2_ y 2 -symmetry of the SC gap 
is the well-established fact [2j, less is known about the 
gap symmetry in HF Cell5s. In the absence of ARPES 
data, the positions of nods in the SC gap for Cell5s are 
subject of debates, and there are only indirect arguments 
in favor of the d x 2_ y i symmetry in CeCoIn^ [121 ]. Fur- 
thermore, CeCoInc, is the multi-band superconductor, in 
distinction to the cuprates, and experiments show that 
some of its Fermi surfaces (FS) are prominently three- 
dimensional (3D) 0, [U . In 3D tetragonal materials, there 
are more choices for the symmetry of the SC gap [l3|. 

The subject of this paper is the analysis of recent neu- 
tron experiments on CeCoIns (Refill) . particularly their 
ability to resolve the issue of the gap symmetry. The ex- 
periments observed the spin resonance in CeCoIn^ at the 
commensurate AFM vector Ql° = (1/2,1/2,1/2), and 
at temperatures below 1.35K (down to the lowest mea- 
sured T = 0A7K). The resonance is very likely related 
to superconductivity, as measurements above T c = 2.3K, 
show only a shallow feature around O.GmeV. 
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The authors of [lj| used the analogy with the cuprates, 
where the spin resonance has been observed [15| at 
Q\P = (1/2, 1/2), and argued that their data can be in- 
terpreted as the evidence that the superconducting gap 
in CeCoInc, is the same d— wave gap with A(p + Qq D ) = 
— A(p), as in the cuprates. This is the case if the res- 
onance is a spin exciton (lt3 |. Below we demonstrate 
that three-dimensionality of Qo, at which the resonance 
has been observed in 14j, makes the analogy with the 
cuprates imprecise. We argue that the resonance un- 
likely comes from 3D FS as we found that a 3D exciton 
is too weak. It may come from Q2D electron 15 band 
in CeCoIns, but then it should be visible not only at 
(1/2,1/2,1/2), but also along the line (1/2,1/2,6) (0 < 
b < 1/2), which apparently is not the case experimen- 
tally. We argue that more likely explanation is that the 
resonance is a magnon-type excitation of /—electrons, 
cleared by a superconductivity. 

To proceed, we first briefly review the situation in the 
2D cuprates. Several candidates for the resonance have 
been proposed, including spin exciton, 7r-resonance, the 
mixture of the exciton and 7r-resonance, and a magnon, 
cleared up by a superconductivity [17] . We focus on the 
spin exciton scenario [3] because of its relation to d x 2_ y 2 
superconductivity. 

The spin-exciton scenario associates the resonance 
with the feedback from the d x i_ y i pairing on the spin 
susceptibility of itinerant fermions near the AFM insta- 
bility. The physics of this effect is not sensitive to the 
details of the model as long as the Fermi liquid concepts 
of the FS and the electron-hole excitations around it are 
preserved. The full dynamic spin susceptibility of itiner- 
ant fermions can be quite generally expressed as 



x(Q,0) = /4 — 



Xo 



r 2 -n(Q,n) 



(i) 



where /is is Bohr magneton, and dimensionless £ is pro- 
portional to the magnetic correlation length (£ _1 = 
signals the onset of AFM order), and H(Q, O) is the po- 
larization operator which has the form 



n(Q,n) = i6g 2 Xo Tj2 J 



d 3 p 
(2tt)3 



(G(p, u)G(p + Q, w + O) + F(p, lu)F+( p + Q,u + O)) 

(2) 

where G, F, and F + are the normal and anomalous 
Gor'kov Green functions, correspondingly, and g is the 
exchange constant for the interaction between conduc- 
tion electrons and their collective spin excitations. Spin- 
exciton scenario assumes that the FS crosses the mag- 
netic Brillouin zone boundary at eight points (hot spots) 
separated by Q = Ql° [the prefactor 16 in ([2]) is 2 * 8, 
where extra 2 is the spin factor] 20]. It further assumes 
that spin excitations with momenta near Qq D are com- 
pletely overdamped in the normal state due to strong 



Landau damping into particle-hole fermionic pairs, such 
that /toiI„((3,S1) « i*yQ, where 7 is proportional to 
9 2 Xo v {Ef) / Ep, where v{Ep) is the density of states 
at the FS (note that g 2 xov(Ep) is dimensionless) The 
real part of n„(Q,0) in the normal state can be safely 
approximated by its static value n n (Q,0) which we ab- 
sorb into £ -2 . In the superconducting state, fermions at 
hot spots become gapped, and the behavior of II(Q, O) 
changes. Expanding near the FS and integrating in Eq. 
([2|) over the transverse momentum component, we ob- 
tain Hi 



n sc (Q,o) = l -£[du> 



1 - 



UJ+LJ- + A1A2 



-A? 



A 2 
^2 



(3) 

where uj± — lo ± 0/2 + i5sign(u> ±0/2), and Ai and A2 
are the gaps at the two hot spots separated by Q. 

The remarkable feature in 2D that distinguishes be- 

pairings comes about 
if this factor is nega- 



tween s— wave, d xy and d 
from the factor A1A2 in |3]): 
tive, which is the case for the d x 2_ y 2 symmetry, but 
not the other two symmetries, then (i) H(Q, 0) = 0, 
i.e., the value of £ is unaffected by superconductivity, 
and (ii) i?eII sc (Q,0) diverges logarithmically at 2Ai, as 
7A1 log Ai/|2Ai — O], and ImU sc (Q, u>) jumps discon- 
tinuously at O = 2Ai from zero to 7r7Ai. At O < 2Ai, 
ImH sc (Q, O) = while i?eiI sc (Q, O) gradually decreases 
and behaves as 7T70 2 /(8A!) at small O. This behavior 
of ILj C (Q,0) guarantees that x(Q, 0) from ([TJ) has a pole 
somewhere below 2Ai (a spin exciton) [2l|. The exciton 
frequency is close to 2Ai at small £, but progressively 
shifts down as £ increases. At large £, the pole is at 
O = (8Ai/(7T7)) 1 / 2 £ _1 . An important fingerprint of the 
2D spin exciton scenario is a negative momentum dis- 
persion of the peak, which originates from the fact that 
the energy of the exciton must vanish at the momentum 
which connects nodal points on the Fermi surface 

For s— wave gap symmetry, ReH sc (Q,0) is negative, 
ReIl sc (Q, 2Aj) does not diverge, and the calculations 
show [19] that ReH sc (Q, O) remains negative for all O 
in which case spin exciton does not emerge. 

Consider now whether one can pass this considera- 
tion to a three-dimensional CeCoIn^. The dHvA ex- 
periments on Cell5 have found several small 3D FS and 
a large 3D hole FS, whose size is large enough such that 
this FS contains hot lines - contours of FS points con- 
nected by a 3D diagonal Qq D . To verify whether such 
FS may be responsible for the spin resonance, we now 
re-evaluate the staggered spin susceptibility, Eq. |T]), 
for a model system with a spherical FS with a diame- 
ter ro > \/3/4. The hot lines on this FS are specified by 
cos(? + -\/2sin(9sin(</> + 7r/4) = 3/(4ro), where 9 and <j> are 
azimuthal and polar angles for a point on a hot line. We 
considered various symmetries of the pairing gap orig- 
inating from different representations of the tetragonal 
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group, -D4 [13J: (i) one-dimensional representations A\ g 
(A fe cx fc 2 + hi), B lg (A fe oc fc 2 - fc 2 ), B 2g (A fe oc k x k y ), 
A 2g (Afe oc k x k y {k\ — ky)) (the A 2g gap would have too 
many nodes), and (ii) two-dimensional representation E g 
(A cx k z (k x +ik y ), A oc k z k x , and A oc k z (k x + k y )). We 
found that the 3D resonance in %(Q,£1) is still possible 
only if the pairing gap retains the same d— wave, fc 2 — ky 
symmetry as in the cuprates, i.e., at the spherical FS 
Afe = A sin 2 9 cos 2cf> (in 3D notations, A oc ¥2,2 + ¥2,-2, 
where Yi tm (9,<p) are spherical harmonics). 

Compared to 2D case, the calculation of n(Q, fi) us- 
ing 3D version of Eq. ([3]) involves a summation along 
the whole hot line. In the normal state, this does not 
lead to a new physics - we still have Landau overdampcd 
excitations with H(Q, Q) = ij^uQ. In a superconducting 
state, however, the integration along hot lines leads to 
three essential differences with the 2D case: 

• 1 . by symmetry, hot lines necessary cross the direc- 
tions = (37r/4)(2n+ 1), n = 0,1,2,3, where the 
d— wave gap vanishes. This implies that ImIl(Q, f2) 
is finite at any non-zero frequency 

• 2. for an arbitrary point fc^ along a hot line, 
the gaps Afe and Afc + Q are not simply related, 
the condition Afe = — Afe + Q is satisfied only for 
special symmetry points along a hot line. As a 
consequence, ReU(Q,0) oc <£dk(Afe + Afc + Q) 2 at 
Q = Qo D is finite, negative, and of the order of 

73L>A 

• 3. the logarithmic divergence of Rell(q, f2) at a 
3D analog of 2Ai is removed due to gap variation 
along a hot line. A simple calculation shows that 
Rell(q, ft) only has a cusp at f2 ~ 1.55 A. 

In FigfT] we plot real and imaginary parts of the po- 
larization operator, and the full Imx(Q,Q), Eq. (fTJ) for 
two different £. We see that the resonance is still there, 
but it is rather broad, and the intensity is not enhanced 
compared to that of the normal state. This is in contrast 
to the sharp peak observed in the neutron experiment. 
Likely then, the observed peak does not come from 3D 
Fermi surfaces. 

Still, the observed resonance could potentially be a 
spin exciton. As we said, dHvA experiments have found, 
that Cell5 has at least one Q2D FS (a FS for the elec- 
tron 15-band). The last FS is a slightly corrugated par- 
allelepiped with the base close to a square. If this FS was 
strictly cylindrical along the z-axis, the hot lines would 
be parallel to z, and for d x 2_ y 2 symmetry of the gap, 
Afc+Q = —Afe for any k— point along a hot line. The inte- 
gration over z then would be harmclcss, and the excitonic 
resonance would be identical to the one in a purely 2D 
system, where it is sharp. There are two requirements: 
one is that the resonance frequency must be smaller than 
twice the gap maximum. The resonance in CeCoIn§ is 






FIG. 1: Panels a) and b) - real and imaginary parts of the 
polarization operator H(Q, O.) for 3D spherical FS with hot 
lines. At large frequencies, JmII(f2) — Q/A, where A is the 
maximum of the d x i_ y i gap. We set ro = 1. Panels c) and 
d) - the imaginary part of the full spin susceptibility, Eq. 
U]) (in arb. units) in the normal state (dashed) and d— wave 
superconducting state (solid) for £ _2 /(7 3D A) = 0.1 (c), and 
£~ 2 /(73dA) = 0.25 (d). Observe that Imx(Q,0) in a d x 2_ y 2 
superconductor does has a "resonance" peak, but its intensity 
is not enhanced compared to the normal state. 



observed at f2 res = 0.6meV, below T c = 2.3K. If we 
use a conservative estimate 2A max ~ AT C , we do find 
that fires is indeed smaller than 2A max . The other re- 
quirement is that such Q2D FS must be large enough 
to contain hot spots in the xy plane. For this, the area 
encircled by the extreme dHvA orbits, S ex t, (marked as 
0:1,2,3 in Refs.0, @]) must exceed a quarter of the area of 



the Brillouin Zone: Sbz = (27r/a) 2 (a 
CeCoIns). One has: 



2e7r\ 
~ch) 



4.61410~ s cm in 



(4) 



where T is the experimental dHvA frequency. According 
to Refs. 

the experimental ratio AS ex t/ Sbz for 
ai orbit is 1.15 > 1. This implies that the area of a Q2D 
FS is indeed large enough to contain hot spots. 

Such explanation of the resonance in CeCoIn^, how- 
ever, disagrees with the experiments in one important 
aspect. If the resonance is an effective 2D exciton, it 
should be present not only for (1/2,1/2,1/2), but also 
for the whole set of momenta Q = (1/2, 1/2, b), where 
< 6 < 1/2. In 14 1 the resonance intensity is clearly 
peaked at b = 1/2. The absence of strong resonances 
at other momenta may be the consequence of the facts 
that the Q2D FS is not a perfect cylinder as evidenced 
by the existence of 3 different dHvA orbits, that A gen- 
erally varies along z axis, and that 3D FS also contribute 

to n(Q,n) 

Still, the constraints on the excitonic description call 
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for another possible explanations for the resonance in 
CeCoIn^. A potential candidate is the "magnon" sce- 
nario, originally suggested for the cuprates [181 ] . In ap- 
plication to CeCoIn 5 , this scenario assumes that this sys- 
tem is close to a AFM state, and contains quasi-localized 
spins of /—electrons coupled to each other and to con- 
duction electrons. To account for the localized spins, one 
should include into the denominator of Eq. ([1]) an ad- 
ditional term (Sl/wo) 2 which in the hypothetical absence 
of the damping would give rise to a magnon mode with 
SI = oJoC -1 . F° r itinerant fermions, ujq is of order Fermi 
energy, and such term is negligible for uo ~ A. For HF 
materials, wo is much smaller, and a magnon mode well 
may have energy comparable to A. 

This mode should in principle be present both above 
and below T c , but in the normal state it is washed out by 
Landau damping. In a superconductor, Landau damping 
is strongly reduced at SI < 2A, and the magnon may 
"clear up" . Furthermore, in 3D, ReIl(Q, SI) is negative at 
Q = Qq D i what effectively reduces £ and hence increases 
the stability of SC state against the onset of AFM. 

Parameterwise, a magnon is overdamped above T c if 
7315 > ~ £~ 2 /A. Using 73D - g 2 XQv{E F )/E Fl 

g 2 XQv{Ep) ~ 1, and taking Ep ~ Tjc for HF materials 
(Tk is a Kondo temperature), we find that the inequality 
amounts to (A/Tr-) > £~ 2 , which can well be satisfied in 
Cell5 due to its proximity to AFM. Below T c , the damp- 
ing is reduced and ReH(Q, SI) appears instead. This i?ell 
scales with frequency as 73£>S1 2 /A, such that the denom- 
inator of Eq. JT]) contains two SI 2 terms, and becomes 
(roughly) 



the points on the FS for which d x 2_ y 2 gap vanishes. For 
such Q, the frequency of the exciton must vanish. On the 
other hand, if the resonance is a magnon, its frequency 
should monotonically increase with deviations from Qq D . 

To summarize, in this paper we addressed the issue 
whether the resonance neutron peak observed in the su- 
perconducting state of CeCoInc, may be a spin exciton 
peculiar to d x 2_ y 2 symmetry of the pairing gap. We ar- 
gued that the peak is unlikely a 3D exciton originating 
from 3D FS of CeCoIn^, as such peak is too weak and 
does not raise above the normal state result. The sharp 
excitonic resonance still might come from Q2D FS of 
CeCoIn^, but such resonance should be observed not 
only at Ql° = (1/2, 1/2, 1/2) but along the whole line of 
(1/2, 1/2, b). We argued that the experimental data 14| 
do not immediately support this scenario. We presented 
another, more plausible explanation of the neutron data, 
namely that the resonance is a magnon-type excitation 
of localized /—electrons, restored and modified by su- 
perconductivity. If the resonance is a magnon, its ob- 
servation is not an argument for d x 2_ y 2 gap symmetry. 
We suggest to study the dispersion of the resonance to 
distinguish between the two scenarios. 

We are thankful to Z. Fisk for the discussion and for 
attracting our attention to Ref. [3] • AVC was supported 
by nsf-dmr 0604406. The work of LPG was supported by 
the NHMFL through NSF Cooperative Agreement No. 
DMR-9527035 and the State of Florida. 



f~ 2 - fl 2 



_ l3D ) 

^ A 



(5) 



where £~ f 2 f = f~ 2 - i?ell(0,0). Using u Q ~ A£ and 
7315 > £~ 2 /A, we see that the two SI 2 contributions are 
of comparable magnitude, i.e., the restored resonance is 
partly a magnon, partly an exciton. We also note that, 
for the model considered above, ReH(Q, SI) changes sign 
and is quite small at SI ~ A (see Fig. lb). If the magnon 
energy is in this range, the feedback electronic contribu- 
tion is negligible, and the resonance below T c is simply a 
magnon, restored by superconductivity. 

Quite obviously, the magnon scenario is not peculiar to 
d x 2_ y 2 gap symmetry, and holds also if the superconduct- 
ing gap has an s— wave symmetry. The only difference 
is that, in the s— wave case, ReH(Q, 0) does not change 
sign for Q = Qq D , and remains negative for all frequen- 
cies 0- 

A possible way to distinguish between magnon and 
spin-exciton scenario is to study the dispersion of the 
resonance peak, S1(Q ^ Qo°)- in the spin- exciton sce- 
nario, the dispersion must be negative, at least in some 
range of Q, because once Q becomes equal to the diame- 
ter of the FS (modulo 2tt), it must, by symmetry, connect 
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